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Artist rendering of MSL entry 
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X-ray micro-tomography 


Berkeley Labs, Flicka 


¢ Advanced Light Source (ALS) at the Lawrence 
Berkeley Natl. Laboratory 

¢ Synchrotron electron accelerator used to 
produce 14KeV X-rays 

¢ Used for many research areas, including | 
optics, chemical reaction dynamics, biological E ar 
imaging, and X-ray micro-tomography. 


X-ray micro-tomography 


Collect X-ray images of the sample as you rotate Use this series of images to 
it through 180° “reconstruct” the 3D object 
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Visualization of FiberForm in PUMA at multiple scales 


Porous Materials Analysis (PUMA) 
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Technical Specifications 


¢ Written in C++ 

¢ GUI built on QT 

¢ Visualization module based on 
OpenGL 

* Parallelized using OpenMP for 
shared memory systems 
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Tortuosity Factors 


Quantifies a materials resistance to a 
diffusive flux 

Important in modeling diffusion/reaction 
systems — such as ablative TPS response 
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y= 


7 = tortuosity factors 

€ = porosity 

Dyer = reference diffusion coefficient 
Der s = effective diffusion coefficient 


Surface rendering of FiberForm tomography in PUMA 
V2.1. Visualization contains = 500 million triangles. 


Knudsen Number 


¢ Non-dimensional number which 
defines the diffusion regime 


K A Mean Free Path 
n= >= —_——————————— 
ly Characteristic Length 


¢ Continuum: Kn << 1 
¢ Transitional: Kn ~ 1 
e Rarified: Kn >>1 
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2D diffusivity simulations using a random walk method in 
PuMA. Particle paths are visualized in red. 


Tortuosity Factors 


Quantifies a materials resistance to a 
diffusive flux 

Important in modeling diffusion/reaction 
systems — such as ablative TPS response 


P Dre f 
Der f 


y= 


7 = tortuosity factors 

€ = porosity 

Dyer = reference diffusion coefficient 
Der s = effective diffusion coefficient 


Surface rendering of FiberForm tomography in PUMA 
V2.1. Visualization contains = 500 million triangles. 
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Reference Diffusion Coefficient 


Dyer = reference diffusion coefficient 


Continuum Free Molecular 


Dre = Dyk Doutk does not exist 


© Dene = - 0A, which is undefined as the 
mean free path approaches infinity 

* Dregs therefore must be based on a length 
scale. In this case, the Diffusion coefficient 
through a capillary of diameter [p 


Vv = mean thermal velocity 
A = mean free path 
NL SS 
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Bosanquet Approximation 


* Used to approximate D,..¢ based on known 
values for D, and Dx. [1] 


1 es 
Dreg Dy Dy 


¢ Rewritten for single species diffusion in a 
capillary, D,-¢¢ becomes [2] 


; 1 /Al 
=> -vypi=- 
oe 32 1 


[1] Tomadakis, 1998 
[2] Pollard, 1948 12 


Choice of Length Scale 


1. Define Ip based on an approximate 
geometric length scale for the 


, ‘ 4€ 
material. Typically = oF mean 


intercept length. (Tomadakis, 
Lachaud, Geodict) 

2. Define lp after the simulations are 
complete as the value which makes 
the tortuosity factor vs. Knudsen 
number plot converge to a single 
value. (Zalc) 


0 60 120 um 


Pore size distribution, computed in GeoDict, of FiberForm. 
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Length Scale Option #1 S 


5 


Define lp) based on an approximate geometric 
length scale for the material. Typically “ or mear 
intercept length. (Tomadakis, Lachaud, Geodict, 


1—D Fibers —— :£q.(9) 
e=0.70 — -:£q.(100) 


PuMA) - 7 
- 3 
Bi perpendicu 
. . ry ir 
¢ Most often used in the literature and S to the Fibers 
software & . 
¢ Requires values of 7,,7; and lp in order to en » Flow Parole 
apply the Bosanquet approximation ; ‘S65 7 to the Fibers 
* 7 is no longer a purely geometrical property, Wik sieunaens 
as it is now a function of the Knudsen numbe! 0 le Pa lr he Pe 
¢ Since 7, had no physical meaning without [p, , KNUDSEN NUMBER. Kn 
this can produce confusing results of 7;,< 1 . | 
Figure from Tomadakis, 1993 
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Length Scale Option #2 


Define |, after the simulations are 
complete as the value which makes the 
tortuosity vs. Knudsen number plot 
converge to a single value. (Zalc, PUMA) 


¢ Requires only one value of 7 and a 
computed length scale, lp, in order 
to apply the Bosanquet 
approximation 


Tortuosity Factor 


F . 0 
* 7 is now longer a purely geometrical = ' 3 = , - = = 
F Knudsen Numb 
property, no longer a function of Kn oe 
. . Tortuosity factor vs Knudsen Number for 1D fibers, 
e 
Easier to understand and implement computed in PUMA, showing the parallel and perpendicular 


tortuosity factors for Option #1 and Option #2 
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Choice of Length Scale 


Define lp) based on an approximate 
geometric length scale for the 


, . 4€ 
material. Typically = oF mean 


intercept length. (Tomadakis, 
Lachaud, Geodict) 

Define lp after the simulations are 
complete as the value which makes 
the tortuosity factor vs. Knudsen 
number plot converge to a single 
value. (Zalc) 


0 60 120 um 


Pore size distribution, computed in GeoDict, of FiberForm. 
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Applying Tortuosity Factors 


* Used to compute Der within a porous 
media, with known tortuosity factor, 7, 
known length scale, lj), and known gas 
properties. 


Dre f 


rey — = 


¢ Using Bosanquet approximation to 
approximate D,e-, the equation becomes 


D Ee | Alp 
= — VY _———_ 
MM 3n- \A+lp 


Surface rendering of FiberForm tomography in PUMA 
V2.1. Visualization contains = 500 million triangles. 
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Numerical Methods S 


Continuum Rarified 
¢ Can be solved using typical ¢ Must be solved using particle 
numerical methods such as finite methods to account for Knudsen 
volume and finite difference effects 
1. Geodict - Explicit Jump Solver 1. PuMA-— Random walk solver 
2. PuMA-— Explicit Jump Solver 2. Geodict — Random walk solver 
3. TauFactor — Finite Volume solver (Knudsen regime) 
3. SPARTA — Direct Simulation Monte 


Carlo 


(A) es. MATH 


Laboratories 2 [VIARKET 


C MATH 
2 MARKET 


Random Walk Solver 


Particle method for solving diffusion 
Velocity and mean path for each particle based on 
exponential distribution 
Diffuse reflections are used for surface collisions 
Symmetric boundary conditions 

effi ~ 2t 


(€7) is the mean square displacement of the particles 
Mean thermal velocity, 7, and mean free path, A, are 


imposed to simulate the desired gas species and 
conditions. 
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Wall Collisions 


Diffuse reflections used for surface collisions 
Collision detection can be based on isosurface or 
cuberille grid 


(a) (b) 


Isosurface (a) and cuberille (b) approximations of a 
cylinder with radius 3 voxels. 
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Wall Collisions 


a 
; - oa 
¢ Diffuse reflections used for surface collisions = 
Sc . . ~ 
¢ Collision detection can be based on isosurface or 3 
cuberille grid o. 
0) 
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og 
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Knudsen Number 
Percent difference (isosurface vs cuberille) vs Knudsen 
number for three different ideal geometries ” 


Comparison to Literature S& 


Test Case #1 24 


Tomadakis 


Tortuosity Factor 
a 
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VW. abs 


0.001 0.01 0.1 1 10 100 1000 10000 
Knudsen Number 


The 5% error is likely due to the limitations of computing in 1993. 


e i 3 

3D Fibers, oie Simulations by Tomadakis were using only 200 particles and 
¢ Intersecting, isotropic likely on a small dataset. The PuMA simulations were run on 
* 0.6 porosity 200,000 particles for a total walk length of 10,000 times the 


domain length - 
NN .°;||\ | -;------EE 


Comparison to Literature S 


st Case #2 Tomadakis Perpendicular 

Tomadakis Perpendicular 
PuMA Perpendicular 
PuMA Parallel 


Te 


Tortuosity Factor 


eet Le, Se Se ee 


0.001 0.01 0.1 1 10 100 1000 10000 
Knudsen Number 


; The 5% error is likely due to the limitations of computing in 1993. 
© 1D Fibers, 512 x 512 x 256 Simulations by Tomadakis were using only 200 particles and 


* Non intersecting likely on a small dataset. The PUMA simulations were run on 
° 0.7 porosity 200,000 particles for a total walk length of 10,000 times the 
domain length 


Bosanquet Analysis 


Test Case #1 
sy 
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Bosanquet Analysis 


Test Case #2 
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e FiberForm, 0.8 mm3 0.001 0.01 0.1 1 10 100 1000 
¢ Transverse isotropic Knudsen Number 


* 0.89 porosity 
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Direct Simulation Monte Carlo 


DSMC is a particle method to simulate 
transitional and rarified flows with high 
fidelity 

Very computationally expensive, preventing 
large or frequent simulations 

DSMC diffusion simulations conducted in 
SPARTA, developed at Sandia National Labs. 
Pressure varied to change the mean free 
path, and therefore the Knudsen number 


Used as a verification case for the random 
walk method 


DSMC simulation of transitional flow over the Space 
Shuttle. Sparta.sandia.gov 
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Direct Simulation Monte Carlo asa 


Test Case #1 


Percent Difference 


0 
0.1 1 10 100 1000 


¢ 3D Fibers, 512° 
¢ Intersecting, isotropic 
* 0.6 porosity 


Knudsen Number 
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Conclusion and Outlook S 


Implemented finite difference and BaNSISCEREGON 

random walk tortuosity factor solvers into Artificial MIgcasiSMoaraGhy 

PuMA V2.1 Geneva I Mieneiiesaies 
Demonstrated the necessity of using an 

isosurface collision detection for Visualization Material Properties Material Response 
complex 3d media, a capability which Maren Porosity Yecocccior say 
currently only exists in PUMA a — 
Verified random walk model for tortuosity OpenGL Transient Heat 
factors against Direct Simulation Monte Rendering 


Carlo (DSMC) simulations. rrr meaner 
Recommend changing current SONG HEN 
definitions of tortuosity factor to restore Diusviy/ Toruosy 
é uIK an nuasen 
the value as a purely geometrical 
property. BlofoniaryVohins 
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